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some special quantum networks for simulating some quantum systems are 
given out. Specially, it is obtained that the quantum network for the time 
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^^ , evolution operator which can simulate, in general, Schrodinger equation, 
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^ ■ Simulating quantum systems has such a meaning that using a specially designed quantum 

^- system, for example quantum computer, which is called a simulated system, to simulate 



another so-called physical quantum system. In 1982, Richard Feynman first proposed that a 



D^. quantum system would be more efficiently simulated by a computer based on the principles 



k J ' *^^ quantum mechanics rather than by one based on the principles of classical mechanics. 

j^ . This is because that the size of the Hilbert space grows exponentially with increase of the 

number of particles. A full quantum simulating demands the exponential resources on a 
classical computer so that it is in general intractable. Since the discovery by Shor of a 
quantum algorithm for factoring in polynomial time, there has been tremendous activity 
in the field of quantum computation including quantum simulating. For example, Lloyd has 
shown how a quantum computer is in fact an efficient quantum simulator, a In addition. 
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some general ideas and schemes to several special quantum systems were proposed and 
discussed, a At present, quantum simulating mainly performs a simulation of the dynamics. 
Obviously, it is very important how to design a quantum network for simulating 
Schrodinger equation. However, it is a difficult task by using of the known method given 
by Barenco et. al. It seems to me that the elementary difficulties are original from the 
facts that Hamiltonian is a summation operator of the kinetic energy and the potential en- 
ergy, and the time evolution operator is a multiple product operator of a series of successive 
the small enough time evolution operators within a given precision. Moreover, the small 
enough time evolution operator is only approximately unitary in practical. Therefore one 
needs an alternative method which can give such a quantum network that it is easily able 
to construct a whole network for the summation and product of the simultaneous and/or 
successive transformations. In my paper, O I proposed an universal quantum network for 
quantum computing U, which can be called quantum CPU and be defined as: 

2*-l 

QiU)= n exp{(f/„„|m)(n|®J^)-Ci}, (1) 

m,n=0 

where Ir and Ia are the identity matrices in the register space and the auxiliary qubit space 
respectively, and C\ = Ir c\ = Ir ® \1)aa{0\. If the graphics rules for the factor with 
form exp{{Umn\'m') {n\) ■ C\} are given out, the picture of quantum network Q{U) can be 
drawn easily. It is the most important that Q{U) has two very useful properties 

Q(f/l + f/2 + ■ ■ ■ + Ur) = Q{U,)Q{U2) ■ ■ ■ Q{Ur), (2) 

Q{UiU2 ■■■Ur) = lR®lA + Cl,ll[ CaQ{U,) I CaC\, (3) 

where Ca is the so-called '^ Connector'^ defined as Ca = Ir® ca = Ir® |0)aa(1|- It is used 
to the preparing transformed state so that this prepared state can be used in the successive 
transformation. Furthermore, note that there are the relations c\ = c\ = 0; caCa + c^a^a = 
I A, thus Ca and Ca can be thought of as the fermionic annihilate and create operators 



respectively in the auxiliary qubit. In order to give out the realization of QCPU for the 
product of transformation in a form of full multiplication, eq.(|3) can be rewritten as 



Q{UiU2---Ur) = {lR, 



input 



C\[\{CaQ{U,)]CaC, 



(4) 



out 



while the initial state is now prepared as {[^{t))) input ® (|^(i)) ® \^) A)out- Therefore, it 
seems to me this new construction of the universal quantum network is scalable easily. In 
terms of it, it is easy to design, assemble and scale the whole quantum networks to simulate 
some quantum systems such as free particle, harmonic oscillator, one particle in a constant 
field. In special, it is obtained in this letter that the quantum network for the time evolution 
operator that can simulate, in general, Schrodinger equation. 

First, a basic skill to simulate quantum systems is to discretize the wavefunction which 
describes the quantum state in a finite space as the following: 



il){x = a, t) — > ip{xm, t) I m 



L/N 



(5) 



When the quantum system is limited within a box, it ought to impose the periodic boundary 
conditions, that is 



'lp(Xm+N,t) =llj(Xm,t). 



Therefore, the wavefunction in time t can be written as a vector 
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(6) 



\'^{t)) = Y,HXm,t)\Xm), 



(7) 



m=0 



where \xm) form a set of basis with the properties of orthogonality and completeness 
{{xm\xn) = ^mn, J2m=o \^m){xrn\ = 1) in A^-dimensioual Hilbert space. For two particles, it 
is extended as 



Afi-l Af2-1 
mi=0 m2=0 



(8) 



Of course, the extension to high dimensional is similar. In fact, the above equation also 
can describe two dimensional case. (Usually taking A^i = N2, this means that the particle 



moves in a square box). A quantum register with k qubits can express a state in Hilbert 

space at most with N = 2^ dimensionaL For two particles in one dimensional, its Hilbert 

space should be Ni x N2 = 2^'^^^'^. Thus, in a classical computer, we need exponentially 

increasing bits to deal with this quantum state. But in a quantum computer, we see that 

ki + ^2 qubits, or two quantum registers with ki and k2 qubits respectively, are suitable 

to this task. In fact, this is just one of reasons why simulating a quantum system can 

be more rapidly done by use of a quantum computer than by use of a classical computer. 

Generally speaking, for 3-dimensional and n particles, we can use 3n quantum registers to 

store the discretizing wavefunction. For the system of identical particles the initial state of 

the quantum computer has to be chosen symmetrically or anti-symmetrically. For quantum 

field theory, its discretized method can be similar to one in lattice gauge theory. 

It is worth emphasing how to realize the fundamental operators such as coordinates and 

momentum is a key point. In coordinate representation, coordinates are directly written 

as a diagonal matrix whose diagonal elements are Xm- But the momentum has a little 

complication because it is a derivative action as the following 

.di){x,t) .i){x + Ax,t)-i){x) 
p^(x, t) = -1^^ = -1 . (9) 

So the momentum operator can be defined by: 

P = -i 2 ( T ) ^i\^rn){x„,+i\ - |a;„)(a;^_i|), (10) 

^ / m=0 

where, in order to make the momentum operator is Hermian, an average of the left and right 
derivative has been taken and the periodic boundary conditions (P) has been used. Thus 

N-l 



P 



nt)) = E 



. l lp{Xm+i,t) -llj{Xm~i,t) 

'' 2 (L/N) 



\Xm). (11) 



m=0 

Further, the kinetic energy operator (natural unit system h = c = 1) is obtained: 



1 /"N^ ^ 



T = — 

8/2 \L 



(12) 



^i\Xm){Xm+2\ + \Xm) {Xm~2\) " 2/7N 
_m=0 

It is clear that the potential operator is a diagonal transformation if one only considers 
the local potential U{x) or the external field Ve{x), that is 



N-1 



U =^U{Xm)\Xm){Xm\, (13) 

N-l 

Ve=^ Ve{Xm)\Xm){Xm\, (14) 

m=0 

while the two-body local interaction U{xi,X2) can be written as 

L^ / J / J '-^ \Xmij Xm2)\Xmi: Xm2/ \Xmi: Xfri2\- V "/ 

mi r?i2 

It is easy to extend to high dimensional and many particles formally. For example, in 
two particles case, the fundamental momentum operators are pi = p^^' ® I2,P2 = h ® p^'^' ■ 
Then, the kinetic energy and the potential energy operators can be constructed in similar 
way. 

Choosing a short time step At, the time evolution operator e~'^* reads 

fi(T) = (1 + iH/\tf/^\ (16) 

where Hamiltonian H = T + V can be given out as above. 

Now, we have the ability to simulate some the quantum systems. First, let's consider a 
simple example - a free particle in one dimensional. Its Hamiltonian is if = p^/2yU. It is 
convenient to simulate it in the momentum representation. So we can use two quantum net- 
works to simulate its dynamics. One is the quantum network Q{F) for the discreting Fourier 
transformation, the other is the quantum network Q(e'^*) for the diagonal transformation. 
Then, we use the connector Ca to combine them in turn 

QFp(e'^*F) = {In)^nvut ® (c\CAQ{e'''')CAQ{F)CAC^^ • (17) 

V / out 

It is easy to realize the quantum network for a diagonal unitary transformation by the 
QCPU. When the diagonal matrix is e^^ */^^, it follows that 

2*^-1 

Q(e'P^*/2M) = J] exp{(e'P^*/2^|x„)(x„| ® I a) ■ C^}. (18) 

m=0 

While in my paper Ql , the quantum network for quantum Fourier transformation on a 
fc— qubit register has been obtained: 
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Qin = n Q[B{n)HMon] = n n (^M[iB{n)HUo\m){n\ ® I^] ■ C^}, (19) 



n=0 m=0 n=0 



r)k 



Xr. 



where B{n)H = Ut,J=oBA^'^^/i'^ ' ^Wj and Mo„ = |0)( 

Now, a free particle in one dimensional can be simulated. But, the above method is 
perhaps not the simplest. It can be seen in the following, by virtue of the kinetic energy 
quantum network proposed by this letter the free particle easier can be simulate more easily. 

Secondly, for a harmonic oscillator in one dimensional, its Hamiltonian is H = X]n^(^ + 
l/2)|n)(n|. Because it is diagonal, its realization of QCPU is easy to get: 

2*=-! 

Qho = n exp{iuj{m + l/2)t{\x^){xm\ ^ Ia) ■ Ca}. (20) 

Its action can simulate dynamics of a harmonic oscillator. Obviously, simulating a harmonic 
oscillator is much simpler than a free particle since an appropriate representation leads 
in the diagonal Hamiltonian is chosen. This implies that the picture theory in quantum 
mechanics can simplify the problem. For example, one would like to simulate one particle 
in one dimensional under a constant field u. Its Hamiltonian is H = p'^/2^ + u. To do this, 
the interaction picture can be used and then Vj = e^^^^ue'^^^^ = u and ipi{t) = e™*-?/). Just 
like done for the free particle, it can be simulated easily. 

The last, in order to simulate a general quantum system, it is necessary to know how 
to simulate Schrdinger equation. For simplicity, let's first consider a local potential or an 
external field. The quantum networks for the kinetic energy and the potential energy can 
be constructed as: 

Q{T)=l[expl-^f^\ i\xm){xm\Eim,m + 2)(g)lA)-C{\ (21) 



exp <(-—( — ) {\xm){xm\E{m,m 
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Q{V) = n exp{{V{xm)\xm){xm\ ® Ia) " c\} , (22) 



m=0 



where E{m,n) is a general exchange transformation defined inO. 

Only consider a small time step, the realization of QCPU for the time evolution operator 
reads 

Q(f/(At)) = Q{lR)QiiAtT)Q{iAtV). (23) 

For a finite time, take the product of all the time evolution operators at small enough time 
steps and obtain finally the whole quantum network for the time evolution operator: 



Q(e'^*) = (/ ^ 



' input 



[t/At] 

Ci I n CAQmAt)) I CaCI, 



(24) 

- out 



By using of it, Schrodinger equation can be simulated in general. 

For two particles, the extension of method is direct, but it is not efficient enough in the 
use of computing resources if one does it directly. Although the quantum network for the 
kinetic energy operator is obtained by the same method, but the quantum network for the 
two body potential needs all the 2^ x 2^ basic elements. Thus, a better method is first 
to reduce Schrodinger equation to the mass center and the relative coordinate system. In 
the mass center system, we need to simulate a free practice, and in the relative coordinate 
system, we need to simulate Schrodinger equation in a single-body potential. They are both 
able to be done just stated above. Of course, the efficiency problem here is said with respect 
to the comparison among the different quantum algorithms. Because, the quantum network 
is built in quantum parallelism, that is, it acts on all the states at the same time. Therefore, 
with respect to classical computing, it must be efficient. In this sense, the above method to 
simulate Schrodinger equation can be extended to higher dimensional and more particles. 

In conclusion, our QCPU and its realizations can be used to various problems of quantum 
simulating. Likewise, it can be used to various problems of quantum algorithms. However, 
it is important to find the special quantum network which can reach at higher or the highest 
efficiency. In order to do this, we need to use the fundamental laws of physics, specially 
the principles and features of quantum mechanics, for example coordinate system choice, 
representation transformation, picture scheme and quantum measurement theory if we have 



thought a quantum computing task to be a physical process. Moreover, to simphfy the 
reahzation of QCPU for all computing steps and find the optimized decomposition, we 
have to use the symmetry properties of every step W as possible. In this letter, such 
some examples have been given out. Some other examples in quantum algorithm also were 
obtained in my paper. This research is on progressing. 
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